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Abstract. The long standing problem of non perturbative renormalization of a gauge field theoretical 
Hamiltonian is addressed and explicitly carried out within an (effective) light-cone Hamiltonian approach 
to QCD. The procedure is in line with the conventional ideas: The Hamiltonian is first regulated by suitable 
cut-off functions, and subsequently renormalized by suitable counter terms to make it cut-off independent. 
Emphasized is the considerable freedom in the cut-off function which eventually can modify the Coulomb 
potential of two charges at sufficiently small distances. The approach provides new physical insight into 
nature of gauge theory and the potential energy of QCD and QED near the origin. The so obtained 
formalism is applied to physical mesons with a different flavor of quark and anti-quark. The excitation 
spectrum of the p-meson with its excellent agreement between theory and experiment is discussed as a 
pedagogical example. 

PACS. ll.10.Ef - 12.38.Aw - 12.38.Lg - 12.39.-x 



1 Introduction 

When starting in 1984 with Discretized Light-Cone Quan- 
tization (DLCQ) [J and with a revival of Dirac's Hamil- 
tonian front form dynamics 0, all challenges of a gauge 
field Hamiltonian theory were essentially open questions, 
particularly the non perturbative bound state problem, 
the many-body aspects, regularization, renormalization, 
confinement, chirality, vacuum structure and condensates, 
just to name a few. The step from the gauge field QCD 
Lagrangian down to a non relativistic Schrodinger equa- 
tion was completely mysterious. Now we know better py. 
We have learned how to partition the problem and how to 
shape our thinking in four major steps: 



£qcd , 
H hC m = M 2 \>P), 
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v 2 m, 
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We have understood, for example, that the chiral phase 
transition, in which the quarks are supposed to get their 
mass, is not the major challenge. The challenge is to un- 
derstand what happens after the phase transition, at zero 
temperature. The challenge is to understand the spectrum 
of physical hadrons and to get the corresponding eigen- 
functions, the light cone wave functions. 

The light-cone wave functions for a hadron with 
mass M encode all possible quark and gluon momentum, 
hclicity and flavor correlations and, in principle, are ob- 
tained by diagonalizing the QCD light-cone Hamiltonian 



LC 



P+P- 



Pl, where P ± 



pO±pz 



(2) 



in a complete basis of Fock states with increasing complex- 
ity. For example, the positive pion has the Fock expansion: 



\%+) 



(n|7r + ) \n) 



(xi,k±i)\ud) 



udg/n 



(xi,]t±i)\udg) 



representing the expansion of the exact QCD eigenstate at 
scale A in terms of non-interacting quarks and gluons. The 
i n particles in a Fock state (n) have longitudinal light-cone 
momentum fractions Xi and relative transverse momenta 
kj_j, with 



p+ 



k° + M 

p0 + pz 



»=1 



5>< 



The form of ^ n /H{xi, k_u) is invariant under longitudinal 
and transverse boosts; i.e., the light-cone wave functions 
expressed in the relative coordinates X\ and fcj_j are inde- 
pendent of the total momentum (P + , PjJ of the hadron. 
The first term in the expansion is referred to as the valence 
Fock state, as it relates to the hadronic description in the 
constituent quark model. The higher terms are related to 
the sea components of the hadronic structure. It has been 
shown that the rest of the light-cone wave function is de- 
termined once the valence Fock state is known , with 
explicit expressions given in 

The key issue is to overcome the problem of any gauge 
theory, that the unregulated theory exposes logarithmic 
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singularities. The problem of regularization and renorma- 
lization has been solved in the perturbative context of 
scattering theory, but not in the non perturbative con- 
text of a Hamiltonian. It is addressed to in the first two 
sections and applied in the remainder of this paper. 



2 Regularization 

Canonical field theory with the conventional QCD La- 
grangian allows to derive the components of the total 
canonical four- momentum P+ Its front form version j2j 
rests on two assumptions, the light cone gauge A + = 
IB] and the suppression of all zero modes [SllZj- The front 
form vacuum is then trivial. 

I find it helpful to discuss the problem in terms of 
DLCQ P]E] ■ In the back of my mind I visualize an explicit 
finite dimensional matrix representation of the Light-Cone 
Hamiltonian as it occurs for finite harmonic resolution. 
Such one is schematically displayed in Fig. 2 of 0. All of 
its matrix elements are finite for any finite x and k± . 

The problem arises for ever increasing harmonic reso- 
lution, on the way to the continuum limit: The numerical 
eigenvalues are numerically unstable and diverge logarith- 
mically 15,16 , contrary to the calculations in 1+1 di- 
mension £Q; see also actual DLCQ calculations in 3+1 by 
Hiller [HJ. 

The reason is inherent to Dirac's relativistic vertex in- 
teraction (ki,hi\V\k2,h2;k3,hs), in which some particle 
'1' is scattered into two particles '2' and '3' with their re- 
spective four-momenta k and helicities h, see Fig. ^ The 
matrix element for bremsstrahlung, for example, is pro- 
portional to k±, (ki,1 \V\k 2 , t; £+ 1) oc |kj_|, see Table 9 
in [3], when the quark maintains its helicity while irradi- 
ating a gluon with four- momentum k 1 ^ = (xP + , kx, k 3 ). 
Singularities arise typically when squares of such matrix 
elements are integrated over all kj_ as in the integrations 
of perturbation theory. 

The singularities are avoided a priori by vertex re- 
gularization, by multiplying each (typically off-diagonal) 
matrix element with a regulating form factor F: 



(ki,h 1 \V\k 2 ,h 2 ;k 3l h 3 ) => 

(fci,/n|y|fc2,^2;fc 3 ,/i3> F(Q) 



(3) 



It took several years to realize that it is the Feynman 
four- momentum transfer across a vertex, Q 2 — — (fci — 
k 2 ) 2 , which governs any effective interaction. The minimal 
requirement for such a form factor is 



The job would be done by a step function, F(Q) — 0(Q 2 — 
A 2 ). The limit A — » suppresses the interaction all to- 
gether, the limit A — > oo restores the interaction and its 
problems. Any finite value of A 2 restricts Q 2 to be finite 
and eliminates the singularities. But the sharp cut-off gen- 
erates problems in an other corner of the theory and F(Q) 
must be an analytic function of Q, as to be seen below. 



Vertex regularization takes thus care of the ultravio- 
let divergences. The (light-cone) infrared singularities are 
taken care of as usual by a kinematical gluon mass. 

As usual, regularization is not unique and many ways 
can do that. Dimensional regularization, for example, is 
not applicable in a matrix approach which is stuck with 
the precisely 3+1 dimensions of the physical world. Vertex 
regularization should be confronted with the Fock space 
regularization of Lepage and Brodsky jBj, see also 
which has blocked the renormalization aspects for many 
years. It also should be confronted with [5| and JUj. After 
applauding the light-cone approach [Jj], Wilson and collab- 
orators ^0] have attempted to base their considerations 
almost entirely on a renormalization group analysis, but 
no concrete technology has emerged thus far. 



3 Renormalization 

The non perturbative renormalization of the Hamiltonian 
was stuck for many years by the fact that the coupling 
constant g and the regulator function F(Q) multiply each 
other in Eq. . It was always clear that one may add non 
local counter terms |10|. but is was not clear how they 
could be constructed. Progress has come from recent work 
on a particular model 1 !) . which did allow to formulate a 
paradigmatic example in modern renormalization theory. 

Here is the general but abstract procedure. 

Suppose to have solved Eq.J2J for a fixed value of the 
7 'bare' parameters in the Lagrangian, for the coupling 
constant g — go and the 6 flavor quark masses to/ = 
(to/)o, and for a fixed value of exterior cut-off scale A = 
Ao- Suppose further that these 7+1 parameters are chosen 
such, that the calculated M 2 agree with the corresponding 
experimental values. Next, suppose to change the cut-off 
by a small amount 5A. Every calculated eigenvalue will 
then change by SM 2 . Renormalization theory is then the 
attempt to reformulate the Hamiltonian, such, that all 
changes 5M 2 vanish identically. 

The fundamental renormalization group equation is 
therefore: 



dMf 



dMf 



g=go,mj=m fo ,A=Ao 



o, 



(5) 



for all eigenstates i. Equivalently one requires that the 
Hamiltonian is stationary with respect to small 5 A: 



LC 







(6) 



Hence forward reference to (go, to/ , Aq), to the 'renorma- 
lization point', will be suppressed. 

The Hamiltonian can be made stationary by making g 
and the to/ functions of A, by introducing physical cou- 
pling constants and masses, g and to/, respectively, which 
themselves are functions of the bare g and m / , and which 
are functionals of the regulator F — F. The variation of 
i?LC reads then 



dH LC = og 

da 



i 



5m f 



S F S J^ = 0, 
SF 
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with the familiar variational derivatives. However, since g 
and rrif are themselves functionals of F, this reduces to 



SH 



LC 



5F- 



5F 







Eq.Q as the fundamental equation of renormalization 
theory is then replaced by 



— dF 

5F = SA— = , 



(7) 



since the variational derivative of the Hamiltionian with 
respect to the regulator is unlikely to vanish. 

It can be solved by counter term technology, as fol- 
lows. A counter term is added to the Hamiltonian, whose 
interaction has exactly the same structure except that the 
regulator F(Q) is replaced by C(Q). This defines 



F(Q,A) = F(Q,A)+C(Q,A) 



(8) 



subject to the constraint that the counter term vanishes 
at the renormalization point, 



C(Q,A) 







(9) 



A=A 



The fundamental equation (JJJ) defines then a differential 
equation 

dC{Q-A) _ dF(Q; A) 

dA dA 1 { J 

which, in its integral form, includes the initial condition 

dF(Q, s 



C(Q,A) = - fds 



ds 



F(Q,Ao)-F(Q,A) . (11) 



The renormalized regulator function, F = F + C ', 

F(Q,A)=F(Q,A ), (12) 

is manifestly independent of A. By construction, the value 
of A is determined by experiment. 

One should emphasize an important point: In deriving 
Eq. I|12l) , use was made of assuming the regulator function 
has well defined derivatives with respect to A, The theta 
function of the sharp cut-off, however, is a distribution 
with only ill defined derivatives. 

This raises an other important point: If F(Q,A) is a 
function of Q/A other than a theta function, one must 
specify how the function approaches the limiting values of 
Eq.|@}. The case of the 'soft' regulator 



F(Q,A) = 



A 2 



A 2 + Q 2 



(13) 



is only a very special example. In a more general approach 
the soft regulator plays the role of a generating function 



F(Q,A) 



N 



n=l 



gn 
dA" 



A 2 



A 2 + Q 2 



(14) 



The partials A n d n /dA n are dimensionless and indepen- 
dent of a change in A. The arbitrarily many coefficients 
si,...,Sjv are renormalization group invariants and, as 
such, subject to be determined by experiment. 



4 The effective (light-cone) Hamiltonian 

In a field theory, one is confronted with a many-body prob- 
lem of the worst kind: Not even the particle number is 
conserved. For to formulate effective Hamiltonians more 
systematically, a novel many-body technique had to be 
developed, the method of iterated resolvents |5lll2| . whose 
details are not important here. 

Important is that the effective light-cone Hamiltonian 
iZ C LC has the same eigenvalue as the full light-cone Hamil- 
tonian Hlc and that it generates the bound state wave 
function of valence quarks by an one-body integral equa- 
tion in (£,k|_): 



M 2 ip hlh2 (x,kj_) 



ml + kl 



x 1 — x 



ip hlfl2 (x,k ± ) 



1 ^ f dx'd z k'^ , , 

7T 2 ££j ^X{1-X)X'(1-X>) 

a{Q q ) a{Qq)" 



3tt 

x F{Q q )F{Qq) 



2Q 2 q 2Q? 



x [u{k x M)ru{k[,h' 2 )} [v(k' 2 , h' 2 ) 7fl v(k 2 , h 2 )] . (15) 

One has achieved step 2 of Eq.JTl): i/cLcl^gg) = M 2 \ty qq ) . 
Here, M 2 is the eigenvalue of the invariant-mass squared. 
The associated eigenfunction ipiHi^ ( x i ki_) is the probabil- 
ity amplitude (x, hi; l — x, —k±, h 2 \ ] Fqq) for finding the 
quark with momentum fraction x, transversal momentum 
and helicity hi, and correspondingly the anti-quark. 
Expressions for the (effective) quark masses mi and rn 2 
and the (effective) coupling function a(Q) are given in 
|12j . Q q and Qq are the Feynman momentum transfers 
of quark and anti-quark, respectively, and u(ki,hi) and 
v{k 2 , h 2 ) are their Dirac spinors in Lepage Brodsky con- 
vention [Sj , given explicitly in [3] . They arrange themselves 
in the Lorenz scalar spinor matrix 

(hi,h 2 \S\h[,h 2 ) = [u(k u hi)^u(k'i,h'i)} 
x [v(k 2 ,h 2 ) 7 f M v(k2 i h 2 )] 

which is a rather complicated (matrix) function of its six 
arguments x, x' , kj_, k^_, as tabulated in [TC^. Finally, the 
form factors F(Q) restrict the range of integration and 
regulate the interaction. Note that the equation is fully 
relativistic and covariant. 

It should be emphasized that Eq. l|15|) is valid only for 
quark and anti-quark having different flavors [5lll2| . The 
additional annihilation term for identical flavors is omit- 
ted. At present, it is investigated by 13 . It should also be 
emphasized that the same structure was obtained with a 
completely different method, with Wegner's Hamiltonian 
flow equations fn] • In is also shown why the concept 
of a 'mean momentum transfer', Q 2 — \ (Qq + Qf) is a 
meaningful simplification. It allows to replace Eci. (|15|l by 



M 2 ^ hlh2 (x,kj_) 



k 2 ill 



2 



1 

3^2 



E 



X 



1-x 

dx'd 2 k^ ^^(a;',ki)s(Q) 



y/x(l - x)x'(l - x') Q 



tp hlh2 (x,k±) 

2 R(Q) 
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The form factors F(Q) have made their way into the regu- 
lator function R(Q) = F 2 (Q) . Krautgartner et al JH] and 
Trittmann et al have shown how to solve numerically 
such an equation with a high precision. But since the nu- 
merical effort is so considerable, it is reasonable to work 
first with (over-)simplified models, as specified next. 

The Singlet-Triplet model. Quarks are at relative 
rest when kj_ = and x = x, with x = m\ / (mi + m 2 ) . 
An inspection of Eq. (33) in ^I] reveals that for very small 
deviations from the equilibrium values, the spinor matrix 
(hi, h^\S\hi, h' 2 ) is proportional to the unit matrix, 

(hi,h2\S\h' 1 h' 2 ) ^ 4miro 2 6 h% ^ . (17) 

For very large deviations, particularly for kj 2 3> kj 2 , holds 

Q 2 ~ki 2 , and (II \S\ TI) ~ 2kf . (18) 

The Singlet- Triplet (ST) model combines these aspects: 

(h 1 ,h 2 \S\h' 1 ,h' 2 ) =S hl:h > t S h2 y 2 (h 1 ,h 2 \S\h 1 ,h 2 ) ,(19) 



(hiM\S\h x ,h 2 ) {^^ + 2, iovhi 



Q 2 



4miffl2 



for hi 



ho 



(20) 



For anti parallel helicities hi — —h 2 (singlets) the model 
interpolates between two extremes: For small momentum 
transfer Q, the '2' in Eci. (|18f) is unimportant and the Cou- 
lomb aspects of the first term prevail. For large Q, the 
Coulomb aspects are unimportant and the hyperfine inter- 
action is dominant. The '2' carries the singlet triplet mass 
difference: Its value is understood by g s (j — (— |)j = g s , 
with the spin-g factor g s — 2. For parallel helicities hi = 
h 2 (triplets) the model reduces to the Coulomb kernel. 
The model over emphasizes many aspects but its simplic- 
ity has proven useful for fast and analytical calculations. 
Most importantly, the model allows to drop the helicity 
summations which technically simplifies the problem enor- 
mously. A more detailed investigation of the spinor matrix 
can be found in |2l)j . 

The model can not be justified in the sense of an ap- 
proximation, but it emphasizes the point that the '2', or 
any other constant in the kernel of an integral equation, 
leads to numerically undefined equations and thus singu- 
larities. Replacing the function a(Q) by the strong cou- 
pling constant a s = g 2 /Air completes the model assump- 
tions. Hence forward, the overline bars for the effective 
quantitites will be suppressed. 



5 The potential energy 

It is possible to subtract a c-number from ff C LC aud to 
define an effective Hamiltonian H c ff implicitly by 

#eLC = ( TO 1 + m 2) 2 + 2 (™1 + m 2) HeS . (21) 

Its eigenvalues have the dimension of an energy 
H eS \ip) = E\<p) , 



achieving this way step 3 of Eq.QJ. Note that mass and 
energy in the front form, on the light cone, are related by 



M 



(mi + m 2 ) + 2 (mi + m 2 ) E , (22) 



and not by M 2 = (mi + m 2 ) 2 + 2 (mi + m 2 ) E + E 2 , as 
usual. Only if the energy is negligible as compared to the 
quark masses, i.e. only if (E/(mi + m 2 )) 2 <C 1, the two 
relations coincide. 

A rather drastic technical simplification is achieved by 
a transformation of the integration variable. One can sub- 
stitute the integration variable x by the integration vari- 
able k z , which, for all practical purposes, can be inter- 
preted [3] as the z-component of a 3-momentum vector 
p = (k z , kj_). For equal masses mi = m 2 = m, the trans- 
formation is, together with its inverse, 



k z z (x) = (m 2 +k ± 2 ) 



y/m 2 + kj 2 + k 2 
(x 



x(l — x) 



(23) 



(24) 



Inserting these substitutions into Eq. (|l(jfl and defining the 
reduced wave function ip by 



ip hl h 2 (x,k±) = 



A(p) = \l 



y / A(k z ,k 1 
y 'x(l — x 



-Vhth^kz,^) , (25) 



(26) 



leads to an integral equation in the components of p, in 
which all reference to light-cone variables has disappeared. 
Using in addition the ST- model of Eq. lff5) l. Ea.lfTB)! trans- 
lates for singlets identically into 



mVp) 



2^ 



- 4 [m 2 + 

d 3 P ' 



( Am 2 



(27) 



R(Q) 



m 



<P(P') 



with a c — ^ct s . The equation for the triplets is obtained by 
dropping the '2'. In the ST-model, the helicity arguments 
in the wave functions can be suppressed. Applying the 
relation between mass and energy, as given in Ea. (|22(l . 
the equation is converted to 



E<p(p) 



2m 



2ir 2 



-<P(P) 

d 3 p' ( Am 2 



(28) 



VMp)W) V Q 2 



R(Q) 

4m 2 



since the reduced mass for mi = m 2 — m is m r = m/2. 

The first term in this equation, p 2 /2m r , coincides with 
the kinetic energy in a conventional non-relativistic Hamil- 
tonian. This is remarkable in view of the fact that no 
approximation to this extent has been made. The fully 
relativistic and covariant light-cone approach has no rela- 
tivistic corrections in the kinetic energy! 
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Since the first term in Eg. (|28fl is a kinetic energy, the 
second must be a potential energy — in a momentum rep- 
resentation, fn principle, it could be Fourier transformed 
with e~ lpr to a configuration space with the variable r. 
But due to the factor A(p)A(p r ) in the kernel, the result- 
ing potential energy would be non-local, see f.e. |2~T]. 

The non-locality of the potential is certainly math- 
ematically exact. But I do not expect this to generate 
aspects of leading importance, and avoid it by the simpli- 
fication A(p) ee 1, both in Eqs.&ft and J2HJ- 

With A(p) — 1, the mean four momentum transfer Q 2 
reduces to the three momentum transfer q 2 = (p — p') 2 . 
In consequence, the kernel of Eq.J2EJ), 



a ( Am 2 



R(q) 



(29) 



depends only on q 
local function, 



4m 2 ' 

p'. Its Fourier transform is a 



V(t) = J d 3 q e-^ r CT(q) , 



(30) 



which plays the role of a conventional potential energy in 
the Fourier transform of Eq. 128(1 , i. e. in 



E V>(r) 



2m r 



V(t) 



(31) 



Here is the Schrodinger equation from Eq. (JJJ ! Despite its 
conventional structure it is a front form equation, designed 
to calculate the light-cone wave function i/j(t) — > <p(p) — > 
ip q q(x,k±). 

I conclude this section with a subtle point, which needs 
clarification in the future. The simplification A(p) = 1 is 
different from a non-relativistic approximation. The ap- 
proach is certainly valid also for relativistic momenta p 2 3> 
m 2 , particularly Eas. l(28jl and l(30|l . The reason is that 
A{p) occurs only under the integral. There, the large mo- 
menta are suppressed by the regulator, anyway. 



6 The renormalized Coulomb potential 

Hence forward, I restrict consideration to the triplet case, 
i.e. to Coulomb kernels like U(q) ~ R(q)/q 2 . The renor- 
malized Coulomb potential is always finite at the origin, as 
opposed to the conventional ^ -singularity. It is instructive 
to verify this explicitly for two regulators: 



U(q) = 



'27T 2 q 2 q 2 +X 2 ■ 



for the soft cut-off, 
A 2 ), for the sharp cut-off. 



(32) 



The Fourier transform according to Eq. (|29|l gives 
V(r)-- 

where Si is the Integral Sine. Asymptotically holds: 
lim V(r) = 



(1 - e" Ar ), for the soft cut-off, 
-^r — Si(Ar), for the sharp cut-off, 



— — , for the soft cut-off, 

— — , for the sharp cut-off. 



(33) 



(34) 



oscillator ~r 2 






i barrier 






f \ 






V ^ 






/ Coulomb ~l/r 








y i 
__ i 

i 

i 

strong coupling 


weak coupling 



(arbitrary units) h/mc 



radius r 



Fig. 2. Schematic behavior of the renormalized Coulomb po- 
tential, see also the discussion in the text. 



Both cut-offs produce the conventional Coulomb poten- 
tial. Near the origin, however, holds: 



lim V(r) = a c \ 
i — »o 



(Ar) 



2 

(Ar 



for the soft cut-off, 



, , for the sharp cut-off. 



The renormalized Coulomb potential is finite but the con- 
stant is cut-off dependent. Even the r-dependence differs: 
The soft cut-off gives a linear and the sharp cut-off a 
quadratic dependence. 

The cut-off dependence near the origin is one of the 
most important aspects of the present work and has a 
deep physical reason to be discussed below. Recalling the 
discussion in Sec.|3|and replacing the soft cut-off in anal- 
ogy to Eq.JTH) with 



R(q) = 



71 — 1 



A 2 



A 2 



(36) 



gives straightforwardly the generalized Coulomb potential 



V(r) 



r 
r 



n—l 

N 

-l+e- Ar 5>„(rA) 



-Ar 



n=0 



(37) 



with so = 1. This result illustrates an other important 
point: The Laguerre polynomials are a complete set of 
functions. The term added to the -1 in Ea. l|37[) is thus 
potentially able to reproduce an arbitrary function of r. 
The description in terms of a generating function, as in 
Eas. l|3b| or l|37|) . is therefore complete. 

The physical picture which develops is illustrated in 
Fig- [3 In the far zone, for sufficiently large r, the potential 
energy coincides with the conventional Coulomb potential 
— — . Since the potential is attractive, it can host bound 
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5 10 15 20 



Fig. 3. The dimensionless Coulomb potential WN(y, 0, 1, 0) 
is plotted versus the radius parameter y — \r for different N, 
i.e. from bottom to top for N = 4, 5, 6, 7, 8. 

states which are probably those realized in weak binding. 
In the near zone, for sufficiently small r, the potential 
behaves like a power series co+cir+C2r 2 which potentially 
can host the bound states of strong coupling, provided the 
actual parameter values allow for that. In the intermediate 
zone, the actual potential must interpolate between these 
two extremes, since Ea. (f3*7|) is an analytic function of r. 
Most likely this is done by developing a barrier of finite 
height, depending on the actual parameter values. The 
onset of the near and intermediate regimes must occur for 
relative distances of the quarks, which are comparable to 
the Compton wave length associated with their reduced 
mass. If the distance is smaller, one expects deviations 
from the classical regime by elementary considerations on 
quantum mechanics, indeed. 

The large number of parameter in Eq. l|37|) can be 
controlled by the following construction: The coefficients 
s n in Eq. (|37[) are expressed in terms of only three param- 
eters a, 6, and c, by 

la b c 

Sn = 1 + 7 TVi + 7 *\i + 7 ■ ( 38 ) 

n\ [n — 1)' (n — 2)! [n — 3)! 

The first few coefficients are then explicitly 

so = 1 , 

S\ = 1 + a , 

s 2 = \ + a + b , (39) 

s 3 = | + f + b + c , 

As a consequence, the dimensionless Coulomb potential, 

W N (y; a, b, c) = ™ = - ( - 1 + e^V s n y n ) , (40) 

which depends on r only through the dimensionless com- 
bination y = Xr, is at most a quadratic function of y, 

W N (y;a,b,c)=a + by + cy 2 , (41) 

in the near zone, and thus independent of N. The remain- 
der starts at most with power y N+1 . A value of a = c = 




Fig. 4. The invariant mass-squares of all available tt + — and 

p + -states are plotted versus a counting index n = 0, 1, 2, 

The straight lines correspond to M% = Mq +n%, with the value 
X = 1.39 GeV 2 , taken from Anisovich et al. gSJ. The filled 
circles correspond to states which have been seen empirically 
|2 1| . the empty ones correspond to the predictions 1251 . — Plot 
courtesy of Shan-Gui Zhou. 

and b = 1 should therefore yield a linear set of functions 
W]y(y; a, b,c) — y in the near zone. As shown in Fig.|3]this 
happens to be true for surprisingly large values of y, i.e. 
not only for y <C 1 . The value of N essentially controls the 
height of the barrier. Similarly, Wjq{y\ 0, 0, 1) = y 2 gener- 
ates a set of functions which are strictly quadratic in the 
near zone. Again, N controls the height of the barrier, as 
to be seen below in Fig. 



7 Determining the parameters by experiment 

The QCD-inspired model developed thus far has a con- 
siderable number of renormalization group invariant pa- 
rameters, which must be determined once and for all by 
experiment. 

In doing this [231 > we have been inspired by the work 
of Anisovich et al. Enumerating the excited states 
of a hadron by a counting index n — 0,1,2,..., these 
authors have found the linear relation M 2 — M 2 + n\ 
for practically all hadrons. As an example, I present in 
Fig. 0] the spectrum of the ir- and the p- meson. 

The linear relation between mass-squared and energy 
on the light cone, Ea. (|22|) . allows then to conclude that the 
potential energy in the near zone must be a pure oscillator, 

V(r) = -a + \f t r 2 , (42) 

at least to first approximation, and that b = in Ea. (|41|) . 
If one addresses to reproduce the spectra of all flavor off- 
diagonal triplet mesons (pseudo- vector mesons), except 
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-2 1 — ■ — ■ — ■ — ■ — 1 — ■ — ■ — ■ — ■ — 1 — ■ — ■ — ■ — ■ — 1 
5 10 15 

r [fm] 

Fig. 5. The continuous lines display the generalized Coulomb 
potential V(r) — Q c AWjv(Ar; a, 0, c) in physical units as func- 
tion of r, for the values iV = 4, 5, 6 from bottom to top. The 
circles indicate the experimental eigenvalues E n for the p + . 
They agree with the calculated eigenvalues for N = 6, shown 
by the horizontal lines. — The dashed line displays the har- 
monic approximation; the horizontal lines on the left indicate 
the oscillator states. — See the discussion in the text. 

the topped ones, one has to determine 6 parameters: The 
2 constants from the oscillator model, ct and ft, and the 4 
effective flavor quark masses m u — rrid, m s , m c , and nit,. 
To determine them, one needs 6 experimental numbers, 
and I take from |24j : 

M„ d - t0 =0.768, M us - to =0.892, M„ e . to =2.010, 
M ud - tl =1.465, M„ s - U =1.680, M uMo =5.325, 

all in GeV. The notation should be self-explanatory. For 
example, M U( ^ tl refers to the first excited state of the p + . 
The so obtained parameter values are: 



m u 


m s 


m c 


m b 


c* 


ft 


0.218 


0.438 


1.749 


5.068 


0.880 


0.0869 


GeV 


GeV 


GeV 


GeV 


GeV 


GeV 3 



The numbers differ slightly from those in [53], due to 
choosing the empirical data set different from Eq. i|43|) , but 
yield about the same overall agreement with all available 
experimental states of pseudo-vector mesons. 

Reverting the argument, one concludes as in [23] that 
the oscillator model in Ea. l|42|l explains quite naturally 
the systematics found by Anisovich et al. |25|. But one 
can do even better. 



8 Relating the oscillator model to QCD 

The oscillator model in Ea. (|42|l is only the harmonic ap- 
proximation to the QCD-inspired, generalized Coulomb 



potential in Ea. H37|) . Their parameters are related obvi- 
ously by 

c t = -a c \a , 6 = 0, f t = 2a c X 3 c . (45) 

One needs more experimental information to pin down the 
value of a, c and N. Choosing A as the QCD scale, i.e. 

A = 200 MeV , (46) 

one can use the expressions for a(Q) in ^2] to calculate 
a s = a(0) from the measured value of the coupling con- 
stant at the Z-mass M z = 91.2 GeV, 

a(M z ) = 0.118 , thus a s = 5(0) = 0.1695 , (47) 

as to be shown in greater detail in [22]. Having fixed a c — 
|a s and A allows to calculate a and c from c t and f t , i.e. 

a = -19.5, c= 24.0. (48) 

We are thus able to draw the generalized Coulomb po- 
tential V(r) — a c XWiy(Xr; a, 0, c) for different N as done 
in Fig. The 'experimental' eigenvalues Eo — E$ for the 
/9-meson, obtained by means of E n = (M% — 4mJ)/(4m u ), 
see Eq. (|22|I . are also inserted, including the empirical lim- 
its of error. The experimental error 8E p ^ ~ ±0.5 GeV 
(thus <5Af Pj3 ~ ±0.1 GeV) is hypothetical, since M Pj3 is 
not confirmed. Taking it for granted, the lowest possible 
value for N is thus 

N = 6 . (49) 

This completes the determination of all parameters. They 
are universal within the model. I thank Harun Omerj2H] 
for giving me the exact eigenvalues prior to publication. 



9 Summary and Conclusions 

This work is an important mile stone on the long way from 
the canonical Lagrangian for quantum chromo dynamics 
down to the composition of physical hadrons in terms of 
their constituting quarks and gluons, by the eigcnfunc- 
tions of a Hamiltonian. 

As part of a on-going effort, a denumerable number of 
simplifying assumptions had to be phrased for getting a 
manageable formalism 5 . Among them is the formulation 
of an effective interaction by the method of iterated resol- 
vents |12| . but the strongest assumption in the present 
work is probably the simplifying Singlet-Triplet model in 
Sec. ^ As long as the assumption are not proven at least 
a posteriori, one must speak of an approach inspired by 
QCD. It is advantageous, however, to have a sufficiently 
simple formalism for penetrating the physical content of 
gauge theory by analytical relations. 

The biggest progress of the present work can be found 
in Sects. [2] and [3| It is related to a consistent regulari- 
zation and renormalization of a gauge theory. The ultra- 
violet divergences in gauge theory are caused less by the 
possibly large momenta of the constituent particles, but 
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by the large momentum transfers in the interaction. In a 
Hamiltonian approach, such as the present, one has not 
much choice else than to chop them off by a regulating 
form factor in the elementary vertex interaction. 

The form factor makes its way into a regulator func- 
tion which suppresses the large momentum transfers in the 
Fourier transform of the Coulomb interaction, see Sec. EJ 
The arbitrariness in chopping off the large momentum 
transfers is reflected in the arbitrariness of the potential 
at small relative distances. It is this arbitrariness which al- 
lows for a pocket in the potential which binds the quarks 
in a hadron. 

The problem is then how to fix this function with its 
many parameters, by experiment. In practice this is less 
difficult than anti-cipated, see Sec. It suffices to deter- 
mine only three parameters, two continuous ones and one 
counting index. 

The potential energy of the present work vanishes at an 
infinite separation of the quarks. This seems be be in con- 
flict with the potential energies of phenomenological mod- 
els |2Z| which rise forever. It also seems to be in conflict 
with lattice gauge calculations [281129) . Is a finite ioniza- 
tion limit in conflict also with 'confinement', i.e. with the 
empirical fact that free quarks have not been observed? — 
The present model prohibits free quarks as a stable so- 
lution, since the sum of the constituent quark masses is 
always larger than the mass of the corresponding hadron 
and a pion. Free constituent quarks would hadronize very 
quickly into bound states. This is different from atomic 
physics with its free constituents, where the binding en- 
ergy is always much smaller than the mass of positronium 
proper. 

The most disturbing aspect of the present work is its 
obvious conflict with lattice gauge calculations 28 ,29 and 
their successes. Several points however should be made: I 
have not checked to which extent a linear term in the 
potential is consistent with the excellent agreement be- 
tween theory and experiment presented in this work. - 
Even with present day computers lattice gauge calcula- 
tions can be extrapolated down to such light systems as 
the 7r or the p only with a head ake. - The calculation of 
the potential energy on the lattice rests on the assump- 
tions of static quarks, of quarks with an infinitely large 
mass. Whether this object is the potential energy to be 
used in a non relativistic Hamiltonian is an open question, 
as well as whether its eigenvalue can simply be added to 
the constituent masses to get the invariant mass of physi- 
cal hadrons. In principle, the relation is justified only only 
for sufficiently small coupling constants. 

The present work opens a broad avenue of further ap- 
plications, among them also the baryons and physical nu- 
clei. But much work must be done in the future before such 
a simple approach as the present must be taken serious. 
It is a first step only. 
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